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We establish a new interstitial-vacancy unbinding transition of the Berezinskii-Kosterlitz- 
Thouless type, transforming the three dimensional pancake vortex lattice of a decoupled layered 
superconductor into a defected solid. This transition is the natural finite-field extension of the 
vortex-anti-vortex unbinding transition establishing the zero-field superfluid stiffness. At finite 
Josephson coupling, the defect unbinding transition turns into a topological decoupling transition. 

PACS numbers; 74.60.Ec, 74.60. Ge 
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The soft vortex matter in high-temperature layered su- 
perconductors exhibits a fascinating rich phase diagram 
with a variety of phase transitions modifying the vari- 
ous intra- and interplanar correlations. These include a 
melting transition from a line solid to a line liquid ^^ , 
a sublimation transition taking the solid into a pancake 
vortex gas |j] , and a decoupling transition ^-0] destroy- 
ing the superconducting coherence between the planes. 
In this letter we concentrate on the decoupled limit where 
the three-dimensional (3D) nature of the vortex system is 
solely due to the long range electromagnetic interaction 
between the pancake vortices. We establish that in the 
crystal phase there is a Berezinskii-Kosterlitz-Thouless 
(BKT) phase transition which separates the pancake- 
vortex lattice from a defected 3D solid with a finite den- 
sity of free interstitials and vacancies, see Fig. 1. With 
finite Josephson interactions this defect-unbinding transi- 
tion also triggers the decoupling of the layers, preempting 
other proposed mechanisms of decoupling 
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FIG. 1. Phase diagram for a weakly-coupled layered super- 
conductor. The solid line marks the defect-unbinding transi- 
tion -Bdcf(r) which transforms into the decoupling transition 
BdcciT) in the case of a finite weak Josephson coupling with 
e < d/X. The dashed line shows the melting line B-^niT) 
calculated within a self-consistent analysis. (We have used 
parameters A(0) ^ 2000 A, d Ri 15 A, and Tc ~ 100 K.) 



The vortex lattice in layered superconductors is "soft" 
and melts at low temperatures; for fields B > B\ = 
$o/A^ (A = penetration depth) the transition is gov- 
erned by the shear interaction between pancake vortices 
in the same layer and thus is close to the two-dimensional 
dislocation-mediated BKT melting in the individual lay- 
ers g . The 3D nature of the transition becomes explicit 
only at low fields B < B\, where the electromagnetic 
(tilt) interaction between pancake vortices in different 
layers dominates over the intralayer (shear) interaction 
due to the planar currents B. Here, we study another 
thermodynamic property of the pancake- vortex lattice - 
its susceptibility to the formation of free defects in the 
form of pancake- vortex interstitials and vacancies. It has 
recently been shown for the pancake-vortex lattice in 
layered decoupled superconductors that the self-energy of 
(and the logarithmic interaction between) vacancies and 
interstitials is strongly screened due to the relaxation of 
the vortex lattice surrounding the defects. We find that 
the corresponding interstitial-vacancy pair creation can 
be understood in terms of a simple vortex-anti-vortex 
pair creation in a superconductor with a vortex-induced 
suppression of the superfluid density. The usual zero-field 
BKT transition establishing the superfiuid stiffness of 
the individual superconducting layers [|lO[ therefore has 
its counterpart at finite fields in the form of interstitial- 
vacancy unbinding in the vortex solid. 

We start from a rigid vortex lattice in a decoupled lay- 
ered superconductor and move one pancake by the dis- 
tance R within the same layer to produce an interstitial- 
vacancy defect pair. The excitation energy associated 
with this manipulation contains a self-, an interaction-, 
and a core-energy part and reads 
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Here, Eq — {^q/AttX}'^ is the basic line energy associated 
with the vortex, d is the layer separation, A denotes the 
planar penetration depth, oq = \/Wq/B is the lattice con- 
stant, and L the macroscopic extent of the layer. The nu- 



mericals rji^jj quantify the core energies of the interstitial- 
and the vacancy-defect and are of order 0.1 |ll||l^ ]. 

Next, we let the lattice relax around the defect. We 
concentrate on a single interstitial and determine the 
screened self-energy - the interaction between the de- 
fects is screened in the same manner. We first ignore 
the electromagnetic coupling to the pancake vortices in 
the other layers and take them into account in a second 
step later. The interaction energy of the pancake-vortex 
configuration reads 
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{Vpc{R) = —2eod\n{R/ao) is the pancake interaction en- 
ergy and Uk = J2i e'^'^' the Fourier transform of the 
pancake- vortex density), from which one can derive the 
dispersive compression modulus 
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with n = B/^o the vortex density. The elastic energy of 
the lattice with one defect can be written in the form 
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where u^ denotes the displacement field of the pancake 
vortex lattice, u^ its Fourier transform. While the first 
term in (^ describes the elastic energy of the relaxed 
lattice, the second term originates from the force the 
interstitial exerts on the pancake vortices, -Esourcc — 
^j y(Ri — VLi) (we place the interstitial at the origin). 
Minimizing (0) with respect to u we obtain the distortion 
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where we have made use of (pi). The displacement 
u(R) = R/27rni?^ removes exactly one pancake vortex 
from the neighborhood of the interstitial and thus pro- 
duces perfect screening B. As a consequence, the log- 
divergent self-energy in ^) is compensated by the gain 
in elastic energy under the relaxation of the lattice: In- 
serting the solution ^) back into (^), we obtain 
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This perfect screening is spoiled by the interaction of 
the pancake vortices in the layer [n — 0) with those in 
the other layers {n ^ 0). Let us view the lattice as be- 
ing made from columns of pancake vortices. The elec- 
tromagnetic interaction of the pancake vortices in the 
n = 0-layer with the rest of the lattice then can be cast 
into the form of a substrate potential 
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where the curvature a^ = (eod/A^) ln(ao/d) follows from 
the electromagnetic energy associated with displacing a 

While the individual 



single vortex in a column |13 
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contribution from each pancake vortex in the column in- 
volves the prefactor d/X and hence is small, the long- 
range nature of the electromagnetic interaction leads to 
a large number (~ A/d) of terms from distant layers. 
Adding the substrate energy V^ub to the elastic energy 
(y) and minimizing, the displacement field (j^ is reduced 
by the factor I + g with Q 
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(i?ci is the lower critical field.) The screening charge is 
reduced to the fraction 1/(1 +5) of a pancake vortex and 
correspondingly the divergent self-energy in (|l|) is only 
partially cancelled by the lattice relaxation; inserting Uk 
back into the expression for the elastic energy, we ob- 
tain Eel — —£od^n{L/aQ)/{l + g). The relaxation of the 
pancake- vortex lattice then reduces the self-energy in (|^) 
by the factor g/{l + g). At large fields, B > iJci, where 
the vortex-vortex interaction is dominated by the intra- 
planar forces we have g < 1 and find a large downward 
renormalization of the energy scale of the defects. When 
5 > 1 at small fields B < Hc-^, the interaction between 
the pancake vortices is dominated by the interlayer forces 
and the effect of screening is small. 

The above analysis extends trivially to the case of a 
pancake-vortex vacancy as well as to the interaction be- 
tween the defects. The energy for an interstitial-vacancy 
pancake-vortex pair of extent R placed into one plane 
of a decoupled layered superconductor then involves the 
screened self- and interaction-energies 
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Before proceeding, we comment on the approximation 
we have made in determining the substrate potential (o), 
where we have ignored the relaxation of the pancake- 
vortex positions in neighboring layers: a simple estimate 
shows that the relaxation Un^i of the pancakes in the 
n-th layer is smaller by the factor d/X, Un,i ~ {d/X)uo^i, 
and hence the correction to the substrate potential (pi) is 
smah, SVsuh,i ^ {X/d)asu'^n,i ^ {d./X)asul^. 

We now discuss the implications of the result (@) for 
the pancake- vortex phase diagram. An immediate conse- 
quence of the logarithmic interaction between the defects 
is a Berezinskii-Kosterlitz-Thouless transition poi. The 
field dependent screening of the interaction pushes this 
defect-unbinding transition to low temperatures. 
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and solving for B we obtain the defect-unbinding line 
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In the presence of a finite Josephson coupling, the appear- 
ance of free defects in the solid triggers the decoupling 
of the superconducting layers and the defect-unbinding 
line B^c[(T) transforms to a decoupling line B^^ciT). For 
fields above B^^dT) the system develops a finite c-axis re- 
sistivity with pc proportional to the number of free mobile 
defects Ud oc a^^ exp(-26/v/l - T/Tdei), see Ref. jl| {b 
is a non- universal constant). This topological decoupling 
transition corresponds to the quartet-unbinding transi- 
tion first proposed by Feigel'man et al. [|l3[ (or the su- 
persolid transitio n pl| ; a defect proliferation transition is 
also described in ]16|]) — here, we establish a rigorous ba- 
sis for the transition via the BKT scenario, accounting for 
the screening response of the pancake-vortex lattice. A 
finite Josephson coupling does not change result (ph if the 
extra Josephson energy ^^gf-Ej within the coherence area 
^dcf ^ ^d remains small, hence a^e^eo/d < Tdcf ~ gsod, 
producing the condition e < d/X. 

In the past, a number of scenarios have been proposed 
for the decoupling transition in layered superconduc- 
tors. We have to distinguish between weak (A < A) and 
intermediate (A < A) coupling, where we have defined 
the Josephson screening length A = d/e. For weak cou- 
pling we define the crossover field Bx — B\h? /}? > B\. 
Below Bx , phase fluctuations are determined by the 
interlayer (tilt) interaction and the decoupling line 
has been predicted to take the form B^cdT) — 
Bx{eod/T:T)l-a{aQ/d) |^,0. Above B^, the intralayer 
(shear) interaction is important and the decoupling line 
takes the form Bdoc(T) ~ B;^{eod/Tf g. For an in- 
termediate coupling the crossover field takes the form 
Bx = Sa = 'I'o/A^ and the low and high field decou- 
pling fines read Bdcc(r) = BAieod/eirT) < i?x ||§1 and 
Sdoc(T) - Ba(£oc?/T)2 > Bx ||. Parametrically, aU 
these results are easily reproduced via a Lindemann cri- 
terion for the interlayer phase correlator MM ■ Our result 
for topological decoupling (^ is valid for weak coupling 
below Bx and so should be compared to the results of 
Daemen et al. [^ and of Horovitz and Goldin [0. Of 
course, whichever decoupling transition occurs first can 
be the only such transition, and we find that our topolog- 
ical decoupling line i?dec(T) is a factor 8 below the line 
advocated in Refs. |6[0|. This numerical relation between 
the two results is not fortuitous, as we will now discuss. 

In fact, the result (S) and the associated transition (0) 
is easily understood as the finite-field extrapolation of the 
zero-field vortex-anti-vortex pair creation and their un- 
binding transition. This understanding brings forward an 
interesting parallel between the competing zero-field in- 
tralayer vortex unbinding- [n3| and the interlayer vortex 
loop transition suggested by Friedel fl^ on the one hand. 



and the present intralayer defect unbinding- and the in- 
terlayer decoupling transition advocated by Daemen et 
al. and by Horovitz and Goldin on the other hand. In or- 
der to understand that we deal with the identical physics 
in both cases we start from the London free energy of a 
layered superconductor, T = J2n ^" + / ^'^'^ B^/8tt with 
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Here, V = {dx,dy), a = 27rA/<I>o is the planar compo- 
nent of the vector potential, (f)n,n+i = fn+i—fn+Jdzaz, 
and e'^ = ra/M < 1 is the mass anisotropy ratio. For 
vanishing interlayer coupling, the individual layers un- 
dergo a zero-field BKT vortex unbinding transition at 
Tbkt — £od/2 [Q. On the other hand, it has been sug- 
gested by Friedel jlj] , that the coupled system might un- 
dergo a decoupling transition triggered by low-energy in- 
terplanar vortex loops. However, as shown by Korshunov 
||l8| , this transition would take place only at higher tem- 
peratures Tl = 4:£od and is therefore preempted by the 
vortex unbinding transition at Tbkt • 

In finite fields the presence of vortices has to be ac- 
counted for. A simple derivation of the appropriate 
free energy starts from the continuous elastic descrip- 
tion of the pancake-vortex system, combining compres- 
sion, shear, and tilt energies. In the incompressible 
limit we can drop the compression term and write the 
phase variable in terms of the displacement field u, 
tpa — (27ri/ao)(uK x 'K)/K'^. The electromagnetic tilt 
energy takes the form geod{V(p)'^ /2Tr with g given by 
(R). Comparing with (lO), the factor g accounts for the 
suppression of the superfluid density due to the presence 
of the vortices. In order to obtain the correct limit for 
_B ^ we have to include the bare superfiuid density. 
This is easily achieved via a calculation of the effec- 
tive penetration depth Acff in the presence of vortices, 
Agjj = A^ -I- B^/Airai,, where ckl denotes the Labusch 
parameter describing vortex "pinning," Vpin,i = ai,uf/2. 
Here, "pinning" is due to the substrate potential (g) and 
hence a,^ — ag/a^d — AnaQg/aQ. Collecting results, we 
obtain the suppression factor for the superfluid density 
Pcff/p = -^Vj^off ~ .9/(1 + .9)- Ii^ the end, the free energy 



replacing (11) in the presence of vortices takes the form 
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For flelds B < Bx = BxK^ /X^ we can drop the shear 
term oc (A(/3)^ in comparison with the electromagnetic 
and Josephson tilt energies. The free energy (n2|) then is 
identical in form to the one studied by Korshunov p8[ ; 
correspondingly, we can discuss two transitions: Turning 
the phase by ±27r we add/remove a vortex from the sys- 
tem and thus recover the energy (ra) for the interstitial- 
vacancy pair creation and the corresponding unbinding 



transition as given by (g). On the other hand, Friedel's 
loop transition takes place at a temperature Tj^(B) = 
4:[g/{l + g)]eod, which is just the result obtained by Dae- 
men et al. (^ and by Horovitz and Goldin [Q . 

The pair of transitions discussed above is not unique 
to the vortex system discussed here — it appears in the 
context of the XY-m.ode\ (with coupling J) subject to a 
symmetry breaking field oc cos[toi^r] (here, ^^ denotes 
the angle of the spin at the position R, m is an inte- 
ger), where the BKT transition at Tbkt = 7rJ/2 com- 
petes with the unlocking or 'roughening' transition at 
Tr = 87rJ/m^ [|9[. A similar scenario shows up in the 
context of the adsorption of a 2D crystal (elasticity C, 
lattice constant a) on a commensurate substrate (lattice 
constant a/m), where the dislocation mediated melting 
transition at T^ = Ca^ /An competes with the depinning 
or 'roughening' transition at Tr = ACa^ jixm^ [^ (in 
both cases the ratio between the two transition tempera- 
tures involves the factor 16/to^, twice larger than in the 
present layered situation). In that sense, Friedel's loop 
transition corresponds to the ijn = 1) roughening tran- 
sition which is preempted by a topological transition at 
lower temperatures. 

Thermal fluctuations will soften the substrate poten- 
tial and modify the result (px|): At high temperatures 
and low fields the renormalized substrate potential (H) 
can be obtained from the thermal average as(T) = 
E„(^^^"("n-"o)/(9uo,x)^)th (seeRef. [|lj for details). 
Going over to Fourier space, the K component of the 
substrate potential is smoothed by the Debye-Waller 
factor exp[— i(r^(u^)th/2]. The mean squared displace- 
ment is in turn determined through the substrate poten- 
tial, (u^)th — IT joLg. Solving self-consistently, we find 
(M2)th = 2AV[(eod/2r)-l)anda,(T) = [£orf/A2(0)](l- 
T/Tbkt)- Here, we have used a mean-field temperature 
dependence \^{T) = A2(0)/(1 - T'^/T^), which gives a 
zero-field BKT transition at Tbkt = £o{T^k.-t) d/2 =4> 
Tbkt = £o(0)d/(H- eo(0)d/Tc). The decoupling line then 
has to be determined from (^ using the renormalized 
suppression factor g{T) = al{l - T/rBKT)/47rA(0)^ for 
the superfluid density and we obtain the final result near 
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The position of the line B'^^^{T) should be compared 
with the melting (or sublimation) line Bn^(T). Using 
the above self-consistent analysis we have determined the 
line in the B-T phase diagram where the substrate po- 
tential collapses and the 3D vortex solid becomes unsta- 
ble pl|. The numerical result is shown as the dashed 
line in Fig. 1, interpolating between Tbkt at zero field 
and the 2D melting temperature T^ « £od/70 at high 
fields. The asymptotic limit at low fields takes the form 
BmiT) = (BA/16^)(rBKT/T)2(l-T/TBKT)2. Comparing 



this with (|l^) we find that the melting line undercuts 
the defect-unbinding line at high temperatures, leaving 
only one (sublimation) transition at B^(T). At low tem- 
peratures, Bdef (T) is below the melting line, allowing for 
two transitions separating a defect-free solid, a phase- 
decoupled solid and a pancake gas. The defect-unbinding 
line B^^f{T) calculated with the self-consistently softened 
substrate is shown as the full line in Fig. 1. 

Though screening strongly reduces the interaction en- 
ergy of the defects for fields B > B\, their core energy 
remains high. For an estimate we consult the numerical 
work of Frey et al. [O and of Olive and Brandt ||l^ , who 
find values of order 0.15 — 0.2eo for the vacancy and inter- 
stitial line defects in a 3D vortex lattice, the interstitial 
typically being 30 % cheaper in energy than the vacancy. 
Therefore at high fields the defect density at the unbind- 
ing transition is small, Ud oc Oq exp[— A^/oq], (with more 
interstitials than vacancies) and the transition is weak. 

In conclusion, we have established a rigorous frame- 
work for the topological decoupling transition in lay- 
ered type II superconductors which is based on a defect- 
unbinding transition of the BKT type. 

We thank the Fonds National Suisse for financial sup- 
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